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This study reports the first measurement of the azimuthal decorrelation between jets with pseu- 
dorapidity separation up to five units. The data were accumulated using the D0 detector during 
the f992-1993 collider run of the Fermilab Tevatron at ^fs = f .8 TeV. These results are compared 
to next-to-leading order (NLO) QCD predictions and to two leading-log approximations (LLA) 
where the leading-log terms are resummed to all orders in a s . The final state jets as predicted 
by NLO QCD show less azimuthal decorrelation than the data. The parton showering LLA Monte 
Carlo HERWIG describes the data well; an analytical LLA prediction based on BFKL resummation 
shows more decorrelation than the data. 
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Correlations between kinematic variables in multijet 
events provide a simple way to study the complex topolo- 
gies that occur when more than two jets are present in 
the final state (|l[ B ||). For example, in dijet events 
the two jets exhibit a high degree of correlation, being 
balanced in transverse energy (Et) and back-to-back in 
azimuth (fa. Deviations from this configuration signal 
the presence of additional radiation. Theoretically this 
radiation is described by higher order corrections to the 
leading order graphs. Using the four momentum transfer 
Q 2 in the hard scattering as the characteristic scale and 
DGLAP (Q) evolution in Q 2 , these corrections have been 
calculated analytically to NLO in perturbative QCD (||; 
H). In addition, they are approximated to all orders by 
using a parton shower approach, like HERWIG ^ for 
example. However, there can be more than one char- 
acteristic scale in the process. Similar to deep inelastic 
lepton-hadron scattering at small Bjorken x and large 
Q 2 , hadron-hadron scattering at large partonic center 
of mass energies (s) may require a different theoretical 
treatment. Instead of just resumming the standard terms 
involving InQ 2 , large terms of the type ln(s/Q 2 ) have to 
be resummed as well using the BFKL technique (||) . Del 
Duca and Schmidt have done this and predict a differ- 
ent pattern of radiation, which results in an additional 
decorrelation in the azimuthal angle between two jets, 
as their distance in pseudorapidity (A77 ~ ln(s/Q 2 )) is 
increased (||). 

In this study, the jets of interest are those most widely 
separated in pseudorapidity (77 = — ln[tan(0/2)], where 
9 is the polar angle with respect to the proton beam). 
The D0 detector (|^) is particularly suited for this mea- 
surement owing to its uniform calorimetric coverage to 
1 77 1 < 4.0. The uranium-liquid argon sampling calorime- 
ter facilitates jet identification with its fine transverse 
segmentation (0.1 x 0.1 in A 77 x A fa. Single particle en- 
ergy resolutions are 15%/ VE and 50% /VE (E in GeV) 
for electrons and pions, respectively, providing good jet 
energy resolution. 

The data for this study, representing an integrated lu- 
minosity of 83 nb _1 , were collected during the 1992-1993 
pp collider run at the Tevatron with a center of mass en- 
ergy of y/s =1.8 TeV. The hardware trigger required a 
single pseudo-projective calorimeter tower (0.2 x 0.2 in 
A7/ x A fa to have more than 7 GeV of transverse energy. 
This trigger was instrumented for \rj\ < 3.2. Events sat- 
isfying this condition were analyzed by an on-line proces- 
sor farm where a fast version of the jet finding algorithm 
searched for jets with Et > 30 GeV. 

Jet reconstruction was performed using an iterative 
fixed cone algorithm. First, the list of calorimeter towers 
with Et > 1 GeV (seed towers) was sorted in descending 
order. Starting with the highest Et seed tower, a preclus- 
ter was formed from all calorimeter towers with 1Z < 0.3, 
where 1Z = y/ A77 2 + Afa was the distance between tower 
centers. If a seed tower was included in a precluster, it 
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FIG. 1. Typical event topology in multijet events. 

was removed from the list. This joining was repeated un- 
til all seed towers become elements of a precluster. After 
calculating the Et weighted center of the precluster, the 
radius of inclusion was increased to 0.7 about this cen- 
ter with all towers in this cone becoming part of the jet. 
A new jet center was calculated using the Et weighted 
tower centers. This process was repeated until the jet 
axis moved less than 0.001 in rj-tj) space between itera- 
tions. The final jet Et was defined as the scalar sum of 
the Et of the towers; its direction was defined using the 
D0 jet algorithm (10), which differs from the Snowmass 
algorithm ([ll]). If any two jets shared more than half of 
the Et of the smaller Et jet, the jets were merged and 
the jet center recalculated. Otherwise, any ambiguities 
in the overlap region were resolved by assigning the en- 
ergy of a given cell in the shared region to the nearest 
jet. Jet reconstruction was over 95% efficient for jets with 
E T > 20 GeV. Jet energy resolution was 10% at 50 GeV 
and jet position resolution was less than 0.03 in both 77 
and <f>. 

Accelerator and instrumental backgrounds were re- 
moved by cuts on the jet shape. The efficiency for these 
cuts was greater than 95%. Based on Monte Carlo sim- 
ulations, residual contamination from backgrounds was 
estimated to be less than 2%. The jet transverse energy 
was corrected for energy scale, out-of-cone showering, 
and underlying event. This correction was based on min- 
imizing the missing transverse energy in direct photon 
events (12). Small pseudorapidity biases (Sij < 0.03), 
caused by the jet algorithm, were also corrected (|l3|). 

A representative multijet event configuration is shown 
in Fig. 0. From the sample of jets with Et > 20 GeV and 
\r)\ < 3.0, the two jets at the extremes of pseudorapidity 
were selected ( J\ and J2 in Fig. |l|) for this analysis. One 
of these two jets was required to be above 50 GeV in Et 
to remove any trigger inefficiency. The pseudorapidity 
difference (A77 = (771 — 772 j ) distribution for events that 
pass the cuts is shown in Fig. |2[ In Fig. ||, the azimuthal 
angular separation, 1 — Afa'ir (A<f> = fa — fa) is plotted 
for unit bins of A77 centered at A77 — 1, 3, and 5. Since 
each distribution is normalized to unity, the decorrelation 
between the two most widely separated jets can be seen 
in either the relative decline near the peak or the relative 
increase in width as A77 increases. 
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FIG. 2. The pseudorapidity interval, A77 = 1 771 — 772 1 , of 
the two jets at the extremes of pseudorapidity. The coverage 
extends to Arj ~ 6. The errors are statistical only. 
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FIG. 3. The azimuthal angle difference, Acj> = tj>i — <j>a, dis- 
tribution of the two jets at the extremes of pseudorapidity 
plotted as 1 - A(/>/tt for A77 = 1, 3, and 5 (0.5 < A77 < 1.5, 
2.5 < Arj < 3.5, and 4.5 < A77 < 5.5). The errors are statis- 
tical only. 



The decorrelation in Fig. |^ can be quantified in terms 
of the average value of cos(7r — A<f>) ([!]). Figure |^ shows 
(cos(7r — A0)) vs. Arj. For the data, the error bars repre- 
sent the statistical and point-to-point uncorrelated sys- 
tematic errors added in quadrature. In addition, the 
band at the bottom of the plot represents the corre- 
lated uncertainties of the energy scale and effects due 
to hadronization and calorimeter resolution. Also shown 
in Fig. I are the predictions from HERWIG, NLO QCD 
as implemented in JETRAD (^|), and the BFKL resum- 
mation (^; The errors shown for the three QCD 

predictions are statistical only. 

The systematic errors, especially the energy scale un- 
certainty, dominate the statistical errors for all Arj except 
for Arj = 5. The jet energy scale uncertainty is estimated 
to be 5%. The resulting uncertainty in (cos(-7r — A</>)) var- 
ied from 0.002 at Ar? = to 0.011 at Arj = 5. Since the 
out-of-cone corrections depended on the pseudorapidity 
of the jet and may not be well understood at large pseu- 
dorapidities, the full size of the out-of-cone showering 
correction was included in the energy scale error band. 
This size of this error was less than 0.013. Uncorrelated 
systematic uncertainties due to the rj bias correction and 
angular resolution were included. This error was less 
than 0.002. The jet selection cuts introduced a system- 
atic uncertainty less than 0.007, which is independent of 
<f> and 77. The uncertainty due to jet position reconstruc- 
tion was estimated by analyzing a subset of the data, 
specifically events with a large Arj, using both Snow- 
mass and D0 jet finding algorithms; the differences in 
(cos(7r — A<f>)) was less than 0.002. 

Comparison of theory with data requires the connec- 
tion of partons with jets. Since no attempt has been 



made to correct the data back to the parton level, the 
the size of the hadronization and calorimeter resolu- 
tion effects were included as an additional systematic er- 
ror. These effects were estimated using HERWIG with 
a detector simulation based on GEANT (H). Jets be- 
fore hadronization were compared with jets after both 
hadronization and detector simulation. In both cases a 
cone jet algorithm with a radius of 0.7 was used. Jets 
reconstructed using partons and particles produced in- 
distinguishable results for (cos(ir — A<fi)) ; the calorimeter 
smearing effects, although negligible for Arj < 3, were 
- 0.02 at A77 = 4 and - 0.03 at Ary = 5. The size of 
these effects were included in the correlated systematic 
error band. 

Since NLO is the first order in perturbative QCD where 
decorrelation is predicted, it may be sensitive to the 
choice of cutoff parameters (scales) necessary in a per- 
turbative calculation. Similar effects have been seen in 
NLO predictions of jet shape (|l6| ) and topologies with 
jets beyond the two body kinematic limit (|l7|). To es- 
timate the size of these effects, the renormalization and 
factorization scales in JETRAD were varied simultane- 
ously from p™ ax /2 to 2p™ ax , where p™ ax is the trans- 
verse momentum of the leading parton. The predic- 
tions for (cos(-7r — Ac/))) varied by less than 0.026. The 
effect of using different parton distribution functions 
(CTEQ2M df), MRSD— (£|), and GRV ©) produced 
variations in JETRAD that were less than 0.0025. Since 
NLO QCD might be sensitive to the jet definition, the 
jet algorithm angle definition study, previously done with 
data, was repeated using JETRAD. The difference be- 
tween the Snowmass and D0 definitions was smaller than 
0.013 for all Arj. 
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FIG. 4. The correlation variable used in this analysis, 
the average value of cos(-7r — A<j>) vs. An, for the data, 
JETRAD, HERWIG, and the BFKL calculations of Del Duca 
and Schmidt. 



The data in Fig. |j show a nearly linear decrease in 
(cos(7r — A<fi)) with pseudorapidity interval. For small 
pseudorapidity intervals both JETRAD and HERWIG de- 
scribe the data reasonably well. JETRAD, which is lead- 
ing order in any decorrelation effects, predicts too lit- 
tle decorrelation at large pseudorapidity intervals. The 
prediction of BFKL leading-log approximation, which 
is valid for large a s Arj, is shown for A77 > 2. As the 
pseudorapidity interval increases, this calculation pre- 
dicts too much decorrelation. Also shown in Fig. [|is the 
HERWIG prediction, where higher order effects are mod- 
eled with a parton shower. These predictions agree with 
the data over the entire pseudorapidity interval range 
(0 < A77 < 5). 

In summary, we have made the first measurement of 
azimuthal decorrelation as a function of pseudorapidity 
separation in dijet systems. These results have been com- 
pared with various QCD predictions. While the JETRAD 
predictions showed too little and the BFKL resummation 
predictions showed too much decorrelation, HERWIG de- 
scribes the data well over the entire Arj range studied. 
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